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A strongly non-integrable system is expected to satisfy the eigenstate thermalization hypothesis, 
which states that the expectation value of an observable in an energy eigenstate is the same as the 
thermal value. This must be revised if the observable is an order parameter for a spontaneously 
broken symmetry, which has multiple thermal values. We propose that in this case the system is 
unstable towards forming nearby eigenstates which yield each of the allowed thermal values. We 
provide strong evidence for this from a numerical study of the 2D transverse-field quantum Ising 
model. 


The eigenstate thermalization hypothesis (ETH) can 
explain how an isolated, quantum many-body system in 
an initial pure state can come to thermal equilibrium 
(as determined by measurements of a specified set of ob¬ 
servables) in finite time JT]-[3] . ETH is expected to hold 
in systems without disorder that are sufficiently far (in 
parameter space) from points of integrability, for observ¬ 
ables that are sufficiently simple (e.g. local) functions of 
the fundamental degrees of freedom. In recent years ETH 
has been the subject of intensive analytic and numerical 
investigations, e.g. [SHU; see |15] for an overview includ¬ 
ing the connection to experimental results in cold atoms 
and other systems. 

The key statement of ETH is that expectation values 
of a relevant observable M in an energy eigenstate |a) 
(of the full many-body hamiltonian H ) take the form 

(a\M\a) = M{E a ), (1) 

where A 4(E) is a smooth function of E and E a is the 
energy eigenvalue. In a system with N 1 degrees of 
freedom, this is enough information to show that A4(E) 
is equal, up to 0(N~ 1 / 2 ) corrections, to the canonical 
thermal average of the operator M, 

TV- Me~ H ! kT , 

M{E)= Tr e~ H / kT [1 + 0(N~^)\, (2) 

where the temperature T is implicitly determined as a 
function of energy E by the usual relation 

Tt He~ H / kT 

E = ^m-■ < 3 > 

A second key statement of ETH is that the off-diagonal 
matrix elements of M in the energy basis, (a\M\j3) with 
a / ft are exponentially small in N. This is needed 
to explain why the diagonal matrix elements of eq. 0 
dominate the instantaneous expectation value of M (in a 
generic time-dependent state) at almost all times, which 
in turn is necessary for thermal equilibrium to be main¬ 
tained once it has been achieved. However this aspect of 
ETH will not be our focus. 

The ETH paradigm must be revisited for a system that 
is capable of exhibiting spontaneous symmetry breaking 


(SSB). Suppose that the observable M is an order param¬ 
eter for a global symmetry. At energies corresponding to 
the broken-symmetry phase, and in the infinite-volume 
limit, we expect the system to have states with the same 
energy but with different values of M (that are related 
by the symmetry). In this case, eq. 0 cannot hold as 
written. We conjecture that, instead, the single smooth 
function JVI(E) is replaced by a multivalued function, 
with one branch for each allowed value of the order pa¬ 
rameter. 

We note that the compatibility of ETH and SSB was 
assumed to hold in [16] , in which the tunnelling dynamics 
of the order parameter were studied in “Schrodinger cat” 
states of a quantum Ising model with disordered infinite- 
range interactions. This model has some special features 
that were the main concern of [16], and so an investiga¬ 
tion of the basic issue in simpler models is warranted. 
Additionally, we focus directly on the equilibrium values 
of the order parameter rather than the quantum dynam¬ 
ics of selected states. 

We therefore turn our attention to a well-known and 
much studied model, the quantum transverse-field Ising 
model with constant nearest-neighbor interactions, spec¬ 
ified by the hamiltonian 

^ = ( 4 ) 

(ij) i 

Here of and af are the usual Pauli matrices at a lattice 
site i, and the first term is a nearest-neighbour sum over 
the links of the lattice. Due to the presence of the trans¬ 
verse field term, this is a fully interacting quantum sys¬ 
tem, and in more than one dimension it is nonintegrable. 
This hamiltonian is invariant under the Zi spin-flip trans¬ 
formation generated by the unitary operator X = EX': 
XHX 1 = H. In two dimensions, in the infinite-volume 
limit, this model exhibits a quantum phase transition at 
a critical coupling g c ~ 3.044 Q3- For g > g c , the ground 
state |0) is unique and satisfies X|0) = |0). The magne¬ 
tization operator 

M = Y,ot (5) 

is odd under the symmetry, A’MZ” 1 = — M, and has 
zero ground-state expectation value, (0|M|0) = 0. For 
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g < g c , the ground state is two-fold degenerate, with 
(0±|M|0±) = ±Mo- The two ground states are related 
by the symmetry, X|0±) = |0 =f). At finite temperature 
with g < g c , there are correspondingly two phases sep¬ 
arated by a second-order phase transition at a critical 
temperature T c . For T < T c , the thermal expectation 
value of M has two values ±A 4(E), where E is related to 
Tbyeq. @. At higher temperatures, A 4(E) vanishes. 
A fluctuation-corrected mean-field computation of Ad (A) 
for infinite volume can be extracted from the results of 

PH- 

At finite volume, more care is required. The energy 
eigenstates |a) are discrete and expected to be nondegen¬ 
erate (for g ^ 0). Each must then also be an eigenstate of 
X with eigenvalue ±1. This implies that (a\M\a) must 
vanish, since M is odd under X. However, for g < g c 
and at energies corresponding to T < T c , we expect the 
energy eigenstates to be unstable to a small symmetry¬ 
breaking perturbation. We therefore modify the hamil- 
tonian by adding M with a coefficient e, 

H -+ H + eM, (6) 

with 1. This explicitly breaks the Zi symmetry by 
an amount that is small compared to the energy scales 
in H. As long as e is not much smaller than the mean 
level spacing (which itself is exponentially small in the 
volume for a large system), we expect the exact energy 
eigenstates to be linear combinations of the unperturbed 
eigenstates with (nearly) equal and opposite expectation 
values of M. In a thermodynamically large system, we 
expect the system to be unstable in this way to an in¬ 
finitesimally small perturbation. 

We investigate the validity of this picture by perform¬ 
ing an exact diagonalization of H on a 4 x 5 lattice with 
periodic boundary conditions for various values of g and 
with e = 10 -3 . This is smaller (by about a factor of 
five) than the mean level spacing across the full spec¬ 
trum, but is still large enough to mix nearby eigenstates. 
In addition to the spin-flip symmetry, there is a dis¬ 
crete translation symmetry (in each cartesian direction) 
and a parity symmetry. We present results on the zero- 
momentum, even parity sector, which has 14,676 states 
(to be compared with 1,048,576 states in the full Hilbert 
space). We find comparable results in other sectors and 
with other (small) values of e. We compute (a\M\a) 
for each state, and compare with the infinite-volume, 
fluctuation-corrected mean-field value AA(E a ) computed 
using the equilibrium methods of jl8j. We use the results 
for coordination number z = 6 with a rescaled value for g 
and M, since the direct results for z = 4 exhibit unphys¬ 
ical features such as a first-order phase transition and 
nonconvexity of the critical curve; this scaling method 
is exact in mean-field theory. This calculation is an es¬ 
sentially uncontrolled approximation, but it provides a 
useful benchmark for our numerical results for individual 
eigenstates at finite volume. 


As a point of comparison, we show results for the g = 0 
model in Fig. |T| In this case, the energy eigenvalues are 
integers ranging from —40 to +32. Each energy eigen¬ 
state is degenerate, and can be chosen to be a simultane¬ 
ous eigenstate of each erf; the magnetization M is then 
obtained by summing these eigenvalues. In this system, 
ETH is clearly not satisfied: for every energy eigenvalue 
below the maximum, there are eigenstates with a range 
of values of M. The large number of magnetization val¬ 
ues for each energy can be understood from the fact that 
states with different net magnetization can result in the 
same nearest-neighbor bond energies, depending on how 
the individual spins are arranged into “droplets” of dif¬ 
ferent sizes. 



FIG. 1. The magnetization as a function of energy for en¬ 
ergy eigenstates of the 2D Ising model on a 4 x 5 lattice with 
periodic boundary conditions and zero transverse field. 

Next we consider a small but finite value for the trans¬ 
verse field coefficient, g = 0.25. Now we find that the 
energy eigenstates are all nondegenerate (even for e = 0), 
as expected. We compute the expectation value of M in 
each eigenstate in the zero-momentum, even-parity sec¬ 
tor, and show the results in Fig. [2] We see that for 
energy greater than roughly zero, the magnetization of 
every eigenstate has been compressed towards a value of 
zero. Below zero energy, the magnetization for each en¬ 
ergy eigenstate has moved closer to one of two possible 
values, one positive and one negative. However, these 
two branches are not sharply defined, indicating that 
ETH is not well satisfied for this value of g. The dashed 
line shows the result of fluctuation-corrected mean-field 
calculation in the infinite volume limit, which is in fair 
agreement with the numerical results. 

Results for g = 0.75 are shown in Fig. [3j There is now 
a clear qualitative difference between the energy eigen¬ 
states in the upper and lower portion of the spectrum. 
Above E ~ —5, all states have a near-vanishing magne¬ 
tization, while below E ~ —22, almost all states possess 
a net magnetization which lies in one of two branches, 
one positive and one negative. Between these two en- 
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FIG. 2. Same as Fig. [l] but with a transverse field with 
coefficient g = 0.25. Only states with zero momentum and 
even parity are shown. Dashed line: fluctuation corrected 
mean-field prediction in the thermodynamic limit, assuming 
ETH. 


FIG. 4. Same as Fig. [3] but with a transverse field with 
coefficient g = 1.5. 
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ergies, the separation between the two branches is less 
pronounced as they merge into a single line at zero mag¬ 
netization. We see good agreement with the fluctuation- 
corrected mean-field calculation. 
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FIG. 3. Same as Fig. [2] but with a transverse field with 
coefficient g = 0.75. 

Results for g = 1.5 are shown in Fig. [4] Two mag¬ 
netization branches can still be discerned, but they are 
much less populated, due to the large transverse field 
lowering the critical energy for spontaneous symmetry 
breaking. Rough agreement with the infinite-volume es¬ 
timate is still seen, though we do not expect our z = 6 
rescaling procedure be as accurate for this larger value of 
9- 

Finally, we show results for g = 3.5 in Fig. [5j Now we 
expect to be in the unbroken phase at all energies, and 
we indeed find zero magnetization for all eigenstates. 

In a much larger system, with values of g large enough 
to be sufficiently far from integrability but less than 
the critical value g Cl we expect that the magnetization 
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FIG. 5. Same as Fig. [3] but with a transverse field with 
coefficient g = 3.5. 


branches would be much better defined, with a spread in 
values that is controlled by the symmetry-breaking coef¬ 
ficient e, as long as e is above a minimum value that is 
exponentially small in the number of sites. 

We also note that at a sufficiently large positive en¬ 
ergy, corresponding to a negative temperature, we expect 
a second phase transition from a disordered phase to an 
antiferromagnetically ordered phase. An order parame¬ 
ter for this phase transition is a staggered magnetization, 
Eq. ([ 5 ]) with a minus sign inserted in the sum on alternat¬ 
ing sites. We have not attempted to study this transition. 

We conclude that our numerical results are fully con¬ 
sistent with the coexistence of spontaneous symmetry 
breaking and eigenstate tliermalization. At system en¬ 
ergies where spontaneous symmetry breaking can occur, 
for an observable M that functions as an order param¬ 
eter, the single smooth function M(E) of eq. |T]) must 
be replaced by a multivalued function. After a small 
symmetry-breaking perturbation is turned on, the ex¬ 
pectation value of M in an individual energy eigenstate 
(a\M\a) will lie on one of these branches, with nearby 
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(in energy) eigenstates yielding expectation values on the 
other branches. 

Systems that exhibit eigenstate thermalization have 
many key physical properties encoded in a single eigen¬ 
state, including some nonlocal properties [2]. Since, 
as we have seen, it is possible to accommodate spon¬ 
taneous symmetry breaking within the ETH paradigm, 
this should extend to critical phenomena at energies near 
a second-order phase transition. Thus it should be pos¬ 
sible, in principle, to extract critical exponents from a 
single eigenstate. Study of this question is currently lim¬ 
ited by small system sizes, but we hope to return to it in 
future work. 
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